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ABSTRACT 

Compressional properties of nuclear matter are studied by using the mean held 
theory with the excluded volume effects of the nucleons. It is found that the 
excluded volume effects make it possible to ht the empirical data of the Coulomb 
coefficient Kc of nucleus incompressibility, even if the volume coefficient K is 
small(~ 150MeV). However, the symmetry properties favor K = 300 ± 50MeV 
as in the cases of the mean held theory of point-like nucleons. 
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One way to determine the incompressibility K of nuclear matter from the 
giant monopole resonance (GMR) data is using the leptodermous expansion[l] of 
nucleus incompressibility K{A,Z) as follows. 

K{A, Z) = K + KsfA-^/^ + + K.Z^A-^!^ + ■ ■ ■ ; / = 1 - 2Z/R, (1) 

where the coefficients and are surface term coefficient, volume- 

symmetry coefficient and Coulomb coefficient, respectively. We have omitted 
higher terms in eq. (1). Although there is uncertainty in the determination of 
these coefficients by using the present data, Pearson [2] pointed out that there 
is a strong correlation among K, Kc and the skewness coefficient, i.e., the third- 
order derivative of nuclear saturation curve. (See table 1.) Similar observations 
are done by Shlomo and Youngblood [3]. 

Table 1 


According to this context, Rudaz et al. [4] studied the relation between 
incompressibility and the skewness coefficient by using the generalized version 
of the relativistic Hartree approximation [5]. The compressional and the surface 
properties are studied by Von-Eiff et al. [ 6 ] [7] [ 8 ] in the framework of the mean 
field approximation of the a-oo-p model with the nonlinear a terms. They found 
that low incompressibility {K k, 200MeV) and a large effective nucleon mass 
Mg at the normal density (0.70 < Mg/M < 0.75) are favorable for the nuclear 
surface properties [ 8 ]. On the other hand, using the same model, Bodmer and 
Price [9] found that the experimental spin-orbit splitting in light nuclei supports 
Mg ^ 0.60M. The result of the generator coordinate calculations for breathing¬ 
mode GMR by Stoitsov, Ring and Sharma [10] suggests K ^ 300MeV. 

In previous papers[ll][12], we have studied the effective nucleon mass Mg, 
incompressibility K and the skewness K' in detail, using the relativistic mean 
field theory with the nonlinear a terms [13] and the one with the nonlinear a and 
uj terms [14]. We found that K = 300 ± 50MeV is favorable to account for K, 
Kc, Kyg and the symmetry energy 04 , simultaneously [11] [12]. It was also found 
that the empirical values of K and Kc in table 1 are not well reproduced by these 
model, if A'/;200MeV [11] [12]. 

On the other hand, the incompressibility K, which is calculated in the frame¬ 
work of the quark-meson coupling (QMG) model is rather small(~ 200MeV) 
[15] [16] [17]. The sub-structure or the finite size effect of nucleons may be impor¬ 
tant in calculating K. In this paper, we study Mg, K, K', Kc and K^s, which can 
be calculated in the framework of nuclear matter with aid of the scaling model 
[ 1 ], by using the relativistic mean field theory with the excluded volume effects 
(EVE) of nucleons [18] [19] [20], and compare the results with the GMR data. We 
also examine whether the QMG result is reproduced by the EVE model. 

We use the relativistic mean field theory with the EVE of nucleons [18] [19] [20]. 
As a Lagrangian, we use the a-oo model with the nonlinear a terms as in ref. [20]. 


2 



( For a while, we restrict our discussions to the symmetric nuclear matter and 
do not consider the p meson effects. ) The Lagrangian density consists of three 
fields, the nucleon tjj, the scalar a-meson 0, and the vector cn-meson V0, i.e., 

Lnctl, = - M)tjj + ^d^4>d>^4> - 

+( 7 , 000 ; F^, = d,V, - d,V„ ( 2 ) 

where m„, ps and are cn-meson mass, a-nucleon coupling and cn-nucleon cou¬ 
pling, respectively. The potential U{(j)) includes a nonlinear cubic-quartic terms 
of the scalar field 0; i.e., 

U (0) = im202 + lfe03 + ic00 (3) 

where rUs is a-meson mass, and b and c are the constant parameters which are 
determined phenomenologically. The Lagrangian density (2) is also the same as 
the one used in the theory of Boguta and Bodmer [13]. 

In the mean field theory of the point-like nucleon, at zero-temperature, the 
baryon density Ppt is given as 


Ppt 


37r^ 


^Fi 


(4) 


where is Fermi momentum and A = 2 in the nuclear matter. In the model 
with the EVE [18] [19] [20], the volume V for the N body system of nucleons in 
configurational space is reduced to the effective one, V — NVn, where 14 is the 
volume of a nucleon. According to this modification for the volume, the baryon 
density p is given by 

= P' 

^ l + Kp'’ 

where p' has the same expression as ppt for the given kp- In a similar way, the 
scalar density is given by 

- ds 
1 + Vnp' 

where p'^ has the same expression as the scalar density of the system of the point¬ 
like nucleons and is given by 


A /- kp + \l k% + 

A = + -)], 


(7) 


where M* is the effective nucleon mass. From the equation of motion for the 
scalar meson, M* is given by 


(pi2 


( 8 ) 
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where Cs = gsM/ms, and ^/gs is the ground-state expectation value of the held 
0. The pressure P and energy density e are also given as 




(9) 

and 

1 + K.p' + 2M2'’ 

(10) 

, respectively, where 

= g^M/m^, 



(11) 

and 

e'{kF,M*) = E*pP -P'. 

(12) 


with Ep = ykp + We remark that the pressure and the energy density of 
free point-like nucleons can be given by eqs. (11) and (12), respectively, if we 
replace M* by the free nucleon mass M. The baryonic chemical potential is also 
given by 

fx = E*p + V^P'ikp, M*) + ^p. (13) 

It is easy to check that the equations (9),(10) and (13) satisfy the thermodynam¬ 
ical identity [18], p = (e -|- P)/p, which yields the following relation, 

M* = ^E% - k% = [{M - ai - VnP'ikpo, M*) - C^po/M^ - ( 14 ) 

at the normal density po, where ai is the binding energy of the normal nu¬ 
clear matter and the subscripts ”0” denotes ”at the normal density”. Since 
P'i^kpo, Mq) > 0, pq > 0 and > 0, it is shown that 

M* < {{M - aif - kl^y/\ (15) 


This condition gives the upper bound for Mq. For example, if we put kpo = 
lAfm~^, we get Mq < 0.94M. (We remark that, due to the modihcation of 
baryon density (see eq. (5)), kpo is larger in the theory with the EVE than in 
the theory of the point-like nucleons, for the given pq. ) 

The incompressibility K at the normal density is dehned as 


^ ^ ^ 9po-^|p=po) 


where e = e/p. In this paper, we also calculate the skewness coefficient, i.e., 
the third-order derivative K' of nuclear saturation curve. In our dehnition, K' is 
dehned as 


, o a e, ( 

K = 3pq—— r|p=pQ — 3po- 


2 IP=P0 


A — 3pq L= 


2 IP=P 0 
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In this definition, large K' means that the equations of state (EOS) is stiff at 
high density. 

There is the eight independent parameters in this model; i.e., M, po, oi, 
Rn{= (3I4/47r)^/^), Os, O^, B and C. In our calculations, we put M = 939MeV, 
Po = 0.17fm“^ and ai = 16MeV. The other five parameters i?„. Os, O^, B and O 
are determined phenomenologically. Besides the two conditions for the saturation 
(i.e., Cq = M — ai and P = 0), if Mg, K and K' are given, we can determine 
the five parameters of the model. Conversely, we can calculate Mg, K and K', 
if the five parameters are given. As is seen in eq. (14), Mg is determined if P„ 
and 0„ are given. Therefore, if we give one (two) quantity (quantities) in P„, 
Ot, and Mg and give two (one) quantities (quantity) in Os, P, O, K and K', the 
other five quantities are automatically determined. We also remark that we put 
P„ = 0 ~ 0.9fm in our calculations. It is difficult to do the calculations at very 
large P„(> 0.9fm), which is close to Po(= {3/(47rpo)}^'^^ ~ l.lfm), i.e., the half 
of the averaged distance between two nucleons in the normal nuclear matter. 

First we give P„, Mg and K, and calculate K'. In fig. 1, we show K' 
as a function of P„ for the fixed values of K and Mg. In the cases of small 
K{< 200MeV), K' decreases in the large P„(^0.75fm) region. We remark that, 
in the case oi K = 200MeI/ and Mg = 0.9M, K' does not decrease much as in 
the other cases in fig. 1(a). K' also decreases in the large Rn region, in the cases 
of large K{> 300MeV) and small Mq(= 0.5M), although the absolute value of 
the decrease is not large as in the cases of small K. In the case of large K and 
the large Mq(= 0.9M), K' hardly decreases. The decrease of K' is related to 
the fact that the coefficient C, which has an important role for determining the 
high density behavior of EOS, becomes (more) negative in the large P„ region. 
C becomes (more) negative (i.e., attractive) to cancel the repulsive effect of the 
EVE, to realize the fixed value of K, as P„ increases. ( The decrease of C 
is conspicuous in the cases of the small Mg, because C must also cancel the 
repulsive effects of the small Mg as well as the effects of the EVE. Also in the 
cases of small K, the decrease of C is conspicuous, because the small K must be 
reproduced. ) As a result, K' becomes smaller in the large P„ region, because 
of the much negative C. It is also seen that K' always increases in the region of 
Pn = 0 ~ 0.7fm, as Rn increases, for the large K{> 300MeV). In these cases, the 
absolute value of the increase is large for the large Mg (= 0.9M). Therefore, K' is 
larger at P„ = O.Sfm, which is often used as the nucleon radius, than at P„ = 0, 
in the cases of the large K and the large Mg. On the contrary, K' is smaller at 
Rn = O.Sfm than at P„ = 0, in the cases of the small K and the small Mg. 

Fig. 1(a),(b) 

After K' is determined, we can also calculate Coulomb coefficient in the 
leptodermous expansion (1), using the scaling model, i.e., using the following 
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equation [1], 


3ql f9K' 


where qei is the electric charge of proton. It is easily seen that Kc becomes more 
negative as K' increases in eq. (18), if K is hxed. Kc is negative and the absolute 
value is large in the EOS, which becomes stiffer at high density. In Eg. 2, we 
show Kc as a function of K with the fixed values of and Mq. In the case 
of Mg/M = 0.5(0.7, 0.9), Kc is smaller (more negative) at = 0.8fm than at 
Rn = 0, if 77/)250(260,160)MeV. The reason is that, in these cases, K' is larger 
at Rn = 0.8fm than at = 0, as is seen in Eg. 1. It is remarkable that, in 
the case of Mg/M = 0.5(0.7, 0.9), the Kc is larger (less negative) at = 0.8fm 
than at = 0, if 77.^250(260,160)MeV. Naturally, the reason is that, in these 
cases, K' is smaller at = 0.8fm than at = 0, as is also seen in Eg. 1. The 
EVE make K' smaller and make Kc less negative for the small K. This change 
of Kc makes it possible to Et the empirical values of Kc for the small K. We also 
remark that the change of Kc by the EVE is opposite to the one by the vector 
meson self-interaction (VSI) [12]. It is reasonable that the repulsive EVE has 
opposite efiects to the ones by the attractive VSI. 


Fig. 2 


If we give 77, Kc and Rn, the other parameters of the model are also deter¬ 
mined. ( We remark that 77' is also uniqnely determined by eq. (18), if 77 and 
Kc are given. ) In table 2, we show the examples of the parameters sets which 
reproduce the empirical values of 77 and Kc in table 1 with several valnes of 7?„. 
For simplicity, we use the average values of 77 and 77^ in table 1. 


Table 2(a),(b),(c),(d) 


We could not End the parameter sets which reprodnce (77, 77c) = (200.0, 2.577)MeV 
at any 7?„, as in the case of 7?„ = 0 [11], where the calculated Kc is always smaller 
than the empirical valne 2.577MeV. This fact is nnderstood as follows. At Rn = 0 
and 77 = 200MeV, the largest Mg has the largest Kc which is closest to the em¬ 
pirical valne. However, as is seen in Eg. 1(a), 77' does not become mnch smaller 
in the cases of 77 = 200MeV and large Mq(^0.9M), even if we use a very large 
7?„(> 0.8)fm. As a result, in those cases, 77c hardly becomes large and does not 
reprodnce the empirical valne. On the other hand, in the cases of 77 = 200MeV 
and Mg ^0.7M, the valne of 77c is mnch smaller at 7?„ = 0 than the empirical 
valne. In those cases, although the EVE makes 77' much smaller and makes 77c 
mnch larger in the large 7?„ region, 77c is still far from 2.577MeV. 

In the case of (77,77c) = (250.0, —0.7065)MeV, the empirical values is repro¬ 
duced by this model in the regions of 7?„ 0 ~ 0.6fm and 7?„ = 0.85 ~ 0.88fm. 


6 



The calculated effective nucleon mass Mg increases as increases in the for¬ 
mer region. The solution of parameter set disappears at = 0.6fm, because 
of the upper bound condition (15) for Mg. The solution appears again at 
Rn = 0.85fm. This reappearance of the solutions is related to the fact that 
K' decreases in the large region. The disappearance of the solution in the 
intermediate region of occurs also in the cases of (i^, Kc) =(300.0, -3.990)MeV 
and (350.0,-7.274)MeV. The solution disappears at Rn = 0.72(0.79)fm and reap¬ 
pears at Rn = 0.86(0.85)fm in the case of {K,Kc) =(300.0, -3.990)MeV ((350.0,- 
7.274)MeV). 

In the cases of (77, 77c) = (150.0,5.861)MeV and (143,3.04)MeV, at Rn = 0, 
there is no parameter set, which reproduce the empirical values. However, we 
could find the parameter sets for {K,Kc) =(150.0,5.861)MeV and for(i7, i7c) = (143, 
3.04)MeV, if we put = 0.81 ~ 0.84fm and = 0.76 ~ 0.87fm, respectively. 
(See EOSIO and EOSll in table 2(d). ) The reason is that, as is seen in fig. 1, 
the EVE makes K' much smaller in the large i?„(^0.75fm) region, in the case of 
K ~ 150MeV. 

We remark that the coefficient C is always negative in the these solutions, as 
is seen in the EOS 10 and EOS 11. C becomes (more) negative (i.e., attractive) to 
cancel the large repulsive effects of the large Rn-, as is mentioned before. Negative 
C may cause the difficulty such as a bifurcation of solution of <l>. However, this 
difficulty is modified by introducing the higher terms of $ which hardly affect 
the nuclear matter properties at the normal density. For example, we add the 
following terms to 77(<h). 


t/56($) 


D 

5M 


$5 ^ 


E 

6ii^ 




(19) 


where D and E are the dimensionless constants. If we put D = —0.00755 and 
E = 0.006, for example, we get the EOS 12 in table 2(d). Although, at p = po, 
the EOS 12 has almost the same properties as in the EOS 11, it has only one 
solution as is seen fig 3, where dU{^)/d^ — ps is shown as a function of $, both 
in the cases of the EOS 11 and of the EOS 12. (We remark dU{^)/d^ — ps = 0 
is a equivalent condition to eq. (8). ) 

We also calculate the volume-symmetry coefficient in the expansion (1). 
Because the p-meson effects are important in the symmetry properties [21] [22], 
we add the standard p-meson terms to the Lagrangian [21] [22] [19]. According to 
the modification, the following term is added to the energy density [21] [22] [19]. 




_ ijp 


^m?p 


{Pp Pn) 


2 _ pRp 


8M2 


Ps: 


Ps Pp Pm 


( 20 ) 


where rrip, Pp, Pp and pn are p meson mass, p-nucleon coupling, proton density 
and neutron density, respectively, and Cp = PpM/rUp. In the theory with the 
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and pn 


( 21 ) 


EVE [19], Pp and pn are given by 


Pp 



1 + VnP' 


Pn 

l + K/’ 


, where the expressions for and p'^ have the same expressions as the densities 
of the point-like proton and neutron, respectively. In the mean held theory, the 
inclusion of the p-meson effects do not affect the saturation conditions and do 
not change the properties such as K and in the symmetric nuclear matter. 
Therefore, the determination of the parameters i?„, Cs, B and C in htting 
the data of K and Kc is not affected by this modihcation of the Lagrangian. 
Using the modihed Lagrangian with the parameter sets in table 2, we calculate 
K^s with aid of the scaling model [1]: i.e.. 


K,.. = K. 


sym 


L 9 


K' 


+ 6 ), 


( 22 ) 


where 


L = 3po 


da4 

dp 


I P=P0 1 


K-sym 


d'^ttA 


dpi- 


2 \p=P0i 


and 


1 (9^e 


P3=0- 


(23) 


The results are also summarized in table 2. In these calculations, we determine 
the p meson coupling pp so as to realize = SO.OMeV at p = po. By comparing 
the table 1 and table 2, it is seen that K = 300 ± SOMeV is favorable to account 
for the empirical values of K, K^. and simultaneously, as in the case of 
/?„ = 0 [11], and as in the case of the mean-held theory with the VSI [12]. This 
unchanged conclusion is related to the fact that K^g is very sensitive to the ratio 
K'/K, which is adjusted to the empirical value in table 2. It seems that this 
feature is not changed drastically, if we use the relativistic mean-held theory and 
the scaling model. 

In the last, we examine whether the QMC result is well reproduced by the EVE 
model. According to ref. [17], K = 200MeV and Mq = 0.906M, if the bag radius 
Rb = O.Sfm is used. Since the value of K' and is not shown in the reference, 
we calculate K' by using the hg. 2 in ref. [17], where the saturation curve of the 
nuclear matter in the QMC model is shown. The result is K' ~ —85MeV. (This 
value corresponds to ~ —3.2MeV, which is somewhat larger than the empirical 
value, 2.577 ± 2.06MeV. ) We search the parameter set of the EVE model, which 
reproduces the K = 200MeV, Mq = 0.906M and K' = — 85MeV. The results are 
shown in table 3. In the parameter set, Rn{= 0.568fm) is somewhat smaller than 
Rb{= 0.8fm). The physical meaning of in the EVE model may be somewhat 
different from that of the bag radius Rb in the QMC model. 


Table 3 



In summary, we have studied the compressional properties of nuclear matter 
by using the relativistic mean held theory with the nonlinear a terms and the 
EVE of the nucleons, under the assumption of the scaling model. We found that 
the EVE yields the possibility to reproduce the empirical values of Kc with the 
small K ~ ISOMeV. However, if we require that Kyg should be reproduced as 
well as Kc, K = 300 ± SOMeV is favorable. It seem to be difficult to change this 
conclusion drastically, in the framework of the relativistic mean held theory and 
the scaling model. 
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Table and Figure Captions 


Table 1 

The sets of the empirical values of K, and K^g. (Shown in MeV.) The sets 
1 ~ 5 is the data from the table 3 in ref. [2]. (According to the conclusion in ref. 
[2], we only show the data in the cases of A' = 150 ~ 350MeV.) The set 6 is one 
example of the data of the table IV in ref. [3]. 

Table 2 

Parameter sets htted for the mean values of the empirical values of K and in 
table 1. The K, K^, K', K^g, L and Kgym are shown in MeV, while are shown 
in fm. In EOS12, D=-0.00755 and E=0.006. (See eq. (19).) 

Table 3 

Parameter sets htted for K = 200MeV, Mq = 0.906M, and K' = —85MeV. 
are shown in fm. 

Fig. 1 iP' as a fnnction of Rn with several values of Mq and K. In the hgure (a) 
((b)), the solid line, the dotted line, and the dash-dotted line are the results with 
Mq/M =0.5, 0.7, and 0.9, respectively, for K = 150(300)MeV. In the hgure (a) 
((b)), the bold solid line, the bold dotted line, and the bold dash-dotted line are 
the resnlts with Mq/M =0.5, 0.7, and 0.9, respectively, for K = 200(400)MeV. 

Fig. 2 iPc as a function of K with hxed values of and Mq. The solid line, 
the dotted line, and the dash-dotted line are the resnlts with Mq/M =0.5, 0.7, 
and 0.9, respectively, in the case of = 0. The bold solid line, the bold dotted 
line, and the bold dash-dotted line are the results with Mg /M =0.5, 0.7, and 0.9, 
respectively, in the case of = 0.8fm. The crosses with error bars are the data 
sets 1 ~ 5 in table 1. The solid squares are the data from the table IV in ref. [3]. 
(For simplicity of the hgnre, we omit the error bars in the latter data. ) 

Fig. 3 dU{^)/d^ — Ps as a fnnction of <h. The solid line and the dotted line 

are the results of EOSll and EOS12, respectively. 
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Set 1 

Set 2 

Set 3 

Set 4 

Set 5 

Set 6 

K 

150.0 

200.0 

250.0 

300.0 

350.0 

143 ± 53 

Kc 

5.861 ±2.06 

2.577 ±2.06 

-0.7065 ± 2.06 

-3.990 ±2.06 

-7.274 ±2.06 

3.04 ±4 

K 

■^^vs 

66.83 ± 101 

-46.94 ± 101 

-160.7 ± 101 

-274.5 ± 101 

-388.3 ± 101 

34 ± 159 


Table 1 


Rn 


Cl 

B 

C 

0.568 

106.17 

16.472 

-2.589x10-^ 

1.049 


Table 3 
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EOS 

1 

2 

3 

K 

250.0 

300.0 

350.0 

Kc 

-0.7065 

-3.990 

-7.274 

K' 

-197 

-94.0 

56.0 

Rn 

0.00 

0.00 

0.00 

MljM 

0.910 

0.831 

0.609 

K 

53.83 

-275.5 

-626.8 

L 

76.61 

78.90 

93.81 

K 

■^^sym 

-29.20 

-24.68 

71.27 


42.392 

144.56 

293.29 

C^v 

18.459 

66.337 

197.85 

B 

-0.5282 

-5.067x10- 

■3 1.595x10-3 

C 

5.071 

0.1028 

-1.686x10-3 


90.31 

84.11 

59.94 



Table 21 

:a) 

EOS 

4 

5 

6 

K 

250.0 

300.0 

350.0 

Kc 

-0.7065 

-3.990 

-7.274 

K' 

-197 

-94.0 

56.0 

Rn 

0.60 

0.72 

0.79 

m;/m 

0.931 

0.919 

0.902 

K 

■^^vs 

77.41 

-244.9 

-600.1 

L 

79.75 

84.81 

91.54 

K 

■^^sym 

-9.193 

24.76 

80.95 


18.121 

55.348 

100.08 


0.42953 

2.0916 

6.7682 

B 

-2.063 

-0.2924 

3.295x10-2 

C 

23.00 

3.728 

0.2168 

r2 

__ 

83.83 

75.39 

66.67 


Table 2(b) 
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EOS 

7 

8 

9 

K 

250.0 

300.0 

350.0 

Kc 

-0.7065 

-3.990 

-7.274 

K' 

-197 

-94.0 

56.0 

Rn 

0.85 

0.86 

0.85 

MljM 

0.573 

0.540 

0.643 

K 

1054 

597.5 

-407.3 

L 

150.5 

164.1 

133.1 

K 

■^^sym 

890.6 

1119 

582.8 


289.99 

302.02 

245.09 

C^v 

181.40 

195.55 

144.14 

B 

2.785x10-3 

2.213x10-3 

3.421x10-3 

C 

-7.121x10-3 

-5.987x10-3 

-9.667x10-3 


11.18 

2.037 

23.87 

Table 2(c) 

EOS 

10 

11 

12 

K 

150.0 

143.0 

143.0 

Kc 

5.861 

3.04 

3.04 

K' 

-260 

-190 

-190 

Rn 

0.81 

0.76 

0.76 

M^/M 

0.574 

0.623 

0.623 

K 

■^^VS 

1988 

1060 

1060 

L 

137.4 

116.9 

116.9 

K 

■^^sym 

668.8 

364.5 

364.5 


307.94 

291.23 

287.83 


189.47 

170.39 

170.39 

B 

3.370x10-3 

4.527x10-3 

3.881x10-3 

C 

-7.261x10-3 

-8.976x10-3 

-5.556x10-3 

r2 

__ 

19.83 

36.22 

36.22 


Table 2(d) 
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